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Abstract: We report on the results from the observations in very high energy band (VHE, _E 7 > 100 
GeV) of the black hole X-ray binary (BHXB) Cygnus X-l. The observations were performed with the 
MAGIC telescope, for a total of 40 hours during 26 nights, spanning the period between June and Novem- 
ber 2006. We report on the results of the searches for steady and variable 7-ray signals, including the first 
experimental evidence for an intense flare, of duration between 1.5 and 24 hours. 



Cygnus X-l is the best established candidate for a 
stellar mass black-hole (BH) and one of the bright- 
est X-ray sources in the sky [6]. Located at a 
distance of 2.2 ± 0.2 kpc, it is composed of a 
21 ± 8 M Q BH orbiting an 09.7 lab companion of 
40 ± 10 M Q [18] in a circular orbit of 5.6 days and 
inclination between 25° and 65° [10]. The X-ray 
source is likely powered by accretion and displays 
the canonical high/soft and low/hard X-ray spectral 
states. The thermal soft component is produced by 
the accretion disk close to the BH, whereas hard X- 
rays are thought to be produced by inverse Comp- 
ton scattering of soft photons by thermal electrons 
in a corona or at the base of a relativistic jet. 

Observations in the soft 7-ray range with COMP- 
TEL [13] and INTEGRAL [7] strongly suggest the 
presence of a higher energy, non-thermal compo- 
nent. Images with the VLBA have shown the pres- 
ence of a highly collimated relativistic jet [16]. 
Cygnus X-l could be a "microblazar", where the 
jet axis is roughly aligned with the line of sight 
[14]. The interaction of the jet with the interstellar 
medium appears to be responsible for a large-scale 
(~ 5 pc diameter), ring-like, radio emitting struc- 
ture [9]. 



Three other binary systems have been detected 
so far in the VHE domain: PSR B1259-63 
[1], LS 1+61 303 [3] and LS 5039 [2]. In 
PSR B1259-63 the TeV emission is thought to 
be produced by the interaction of the relativis- 
tic wind from a young pulsar with the outflow of 
the companion star. Recent results suggest that 
LS 1 +61 303 also contains a non-accreting neu- 
tron star [8], while the situation is not yet clear for 
of LS 5039 [12, 15]. To date, however, there is no 
experimental evidence of VHE emission from any 
galactic BHXB system. 

Cygnus X-l was observed with MAGIC [11] be- 
tween June and November 2006 [4] . The data set 
comprises 40.0 good observation hours from 26 
different nights (see Table 1 for details). A search 
for steady 7-ray signals was performed for the en- 
tire recorded data sample, yielding no significant 
excess. This allows us to establish the first upper 
limits to the VHE 7-ray steady flux of Cygnus X- 
1 in the range between 150 GeV and 3 TeV (see 

Figure 1), of the order of < 1 5% of the Crab 

nebula flux. Given the timescale of the variability 
of Cygnus X- 1 in other energy bands, 7-ray signals 
were also searched for on a daily basis. The results 
are shown in Table 1 . We obtain results compat- 
ible with background fluctuations at 99% CL for 
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Table 1 : From left to right: Modified Julian Date 
of the beginning of the observation, total effective 
observation time (EOT), number of excess events, 
statistical significance of the excess, equivalent 
(post-trial) significance for 26 independent sam- 
ples. A cut SIZE>200 photo-electrons (E 7 > 150 
GeV) has been applied. 
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Figure 1: Differential energy spectrum from 
Cygnus X-1 corresponding to 78.9 minutes EOT 
between MJD 54002.928 and 54002.987 (2006- 
09-24). Also shown are the Crab nebula spectrum 
and the best fit of a power-law to the data and the 
95% CL upper limits to the steady 7-ray flux. 
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Figure 2: Distribution of 6 2 values for the signal 
(dots) and background (histogram) events for an 
energy threshold of 150 GeV. 



all searched samples except for MJD=54002.875 
(2006-09-24). We derive upper limits to the inte- 
gral flux above 150 GeV between 2 and 25% of 
the Crab nebula flux (depending basically on the 
observation time) for all samples compatible with 
background fluctuations. The data from 2006-09- 
24 was further subdivided into two halves to search 
for rapidly varying signals, obtaining 0.5(7 and 
4.9a results for the first (75.5 minutes EOT start- 
ing at MJD 54002.875) and second (78.9 minutes 
EOT starting at MJD 54002.928) samples, respec- 
tively. The post-trial probability is conservatively 
estimated by assuming 52 trials (2 per observation 
night) and corresponds to a significance of 4.1cr. 
The sample corresponding to MJD 54002.928 was 
further subdivided into halves, obtaining 3.2a and 



3.5cr excesses in each. At this point we stopped the 
data splitting process. 

The distribution of 92 for signal and background 
events corresponding to the 78.9 minutes EOT 
sample starting at MJD 54002.928 is shown in 
Figure 2. The excess is consistent with a point 
source located at the position of Cygnus X- 1 . The 
map of excess events around the source is shown 
in Figure 3. A Gaussian fit yields the location: 
a = 19 h 58 m 17 s , S = 35°12'8" with statistical 
and systematic uncertainties of 1.5' and 2', respec- 
tively, compatible within errors with the position 
of Cygnus X-1 and excluding the radio nebula at a 
distance of ~ 8'. The energy spectrum is shown in 
Figure 1. It is well fitted (x2/n.d.f = 0.5) by the 
following power law: dN/(dA dt dE) = (2.3 ± 
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Figure 3: Gaussian-smoothed (a = 4') map of 7- 
ray excess events (background subtracted) above 
150 GeV around Cygnus X-1 corresponding to 
78.9 minutes EOT between MJD 54002.928 and 
54002.987 (2006-09-24). The black cross shows 
the best-fit position of the 7-ray source. The po- 
sition of the X-ray source and radio emitting ring 
are marked by the green star and contour, respec- 
tively. The purple dots mark the directions tracked 
during the observations. Note that the bin contents 
are correlated due to the smoothing. 
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0.6) x 10- 12 (£/1 TeV)- 3 - 2±0 - 6 cm- 2 s- 1 TeV" 1 
where the quoted errors are statistical only. We es- 
timate the systematic uncertainty to be 35% on the 
overall flux normalization and 0.2 in the determi- 
nation of the spectral index. 

The excess from the direction of Cygnus X- 
1 occurred simultaneously with a hard X-ray 
flare detected by INTEGRAL [17], Swift/BAT and 
RXTE/ASM. Figure 4 shows the correlation be- 
tween MAGIC, Swift/BAT and RXTE/ASM light- 
curves. The TeV excess was observed on the ris- 
ing side of the first hard X-ray peak, 1-2 hours be- 
fore its maximum, while there is no clear change 
in soft X-rays. Additionally, the MAGIC non- 
detection during the following night (yielding a 
95% CL upper limit corresponding to a flux ~5 
times lower than the one observed in the second 
half of 2006-09-24) occurred during the decay of 
the second hard X-ray peak. A possible explana- 
tion is that, during the night of 2006-09-24, soft 
and hard X-rays were produced in different re- 
gions. Furthermore, hard X-rays and VHE 7-rays 
could be produced in regions linked by the colli- 



Figure 4: From top to bottom: MAGIC, Swift/BAT 
and RXTE/ASM measured fluxes from Cygnus X- 
1 as a function of time. The left panels show the 
whole time spanned by MAGIC observations. The 
vertical, dotted blue lines delimit the range zoomed 
in the right panels. The vertical red line marks the 
time of the MAGIC signal. 

mated jet, e.g. the X-rays at the jet base and 7- 
rays at an interaction region between the jet and the 
stellar wind. These processes would have different 
physical timescales, thus producing a shift in time 
between the TeV and X-ray peaks. Note that the 
distance from the compact object to the TeV pro- 
duction region is constrained below 2' by MAGIC 
observations and therefore it is unrelated with the 
nearby radio emitting ring-like structure [9]. The 
observed TeV excess took place at phase 0.91 
(phase 0.0 is when the BH is behind the O star). 
Currently, MAGIC observations are available only 
for the night 2006-09-24, which precludes any pos- 
sible analysis of periodicity for the TeV emission. 
The jet scenario, however, has some constraints. If 
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the TeV emission were produced in the jet, well 
within the binary system, the photon-photon con- 
version in the stellar radiation field would be very 
substantial, renderiing unlikely a TeV detection[5]. 
Admittedly, the inclination of the orbit and the an- 
gle of propagation to the observer can affect this. 
Even without an explanation for a TeV flare, it is 
possible that the emission could have originated far 
from the compact object. Interactions of the jet 
with the stellar wind may lead to such a situation. 

In summary, for the first time we have found ex- 
perimental evidence of VHE emission produced by 
a Galactic stellar-mass BH. It is also the first evi- 
dence of VHE gamma-rays produced at an accret- 
ing binary system. Our results show that a possible 
steady VHE flux is below the present IACT's sen- 
sitivity and tight upper limits have been derived. 
On the other hand, we find evidence for an intense 
flaring episode during the inferior conjunction of 
the optical star, on a timescale shorter than 1 day 
with a rise time of about 1 hour, correlated with a 
hard X-ray flare observed by Swift and INTEGRAL. 
These results imply the existence of a whole new 
phenomenology in the young field of VHE astro- 
physics of binary systems to be explored by present 
and future IACT's. 

We thank the IAC for the excellent working condi- 
tions at the Observatorio del Roque de los Mucha- 
chos in La Palma. 
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